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Abstract. We review our analysis of the clustering properties of "Lyman-break" galax- 
ies (LBGs) at redsliift z ~ 3, previously discussed in Wechsler et al. [1]. We examine 
the likelihood of spikes found by Steidel et al. [2] in the redshift distribution of LBGs, 
within a suite of models for the evolution of structure in the Universe. Using high- 
resolution dissipationless N-body simulations, we analyze deep pencil-beam surveys 
from these models in the same way that they are actually observed, identifying LBGs 
with the most massive dark matter halos. We find that all the models (with SCDM 
as a marginal exception) have a substantial probability of producing spikes similar to 
those observed, because the massive halos are much more clumped than the underly- 
ing matter - i.e., they are biased. Therefore, the likelihood of such a spike is not a 
good discriminator among these models. The LBG correlation functions are less steep 
than galaxies today (7 ~ 1.4), but show similar or slightly longer correlation lengths. 
We have extcned this analysis and include a preliminary comparison to the new data 
presented in Adelberger et al. [3]. We also discuss work in progress, in which we use 
semi-analytic models to identify Lyman-break galaxies within dark-matter halos. 



INTRODUCTION 



Recently, Steidel et al. (S98) [2] discovered that the redshift distribution of 
"Lyman-break" galaxies (LBGs) in a pencil beam reveals a large "spike" in the 
LBG distribution near z ~ 3.1. This spike corresponds to a fractional overdensity 
of LBGs of a few hundred percent over a comoving scale of order ~ 10 — 20/i -1 Mpc. 
Since then, they have compiled a sample of more than 600 galaxies in six fields, 
with measured redshifts between about 2.5 and 3.5 [3], [4]. 

At a first glance, such a high peak seems surprising, since it suggests substantial 
nonlinear clustering on rather large scales. In fact, we [1] and other authors [5], [6], 
[7] find from simulations that these spikes arise naturally in a variety of cosmologies, 
provided that there is substantial galaxy biasing. 



Here we compare the clustering properties of the LBGs to those expected in high- 
resolution simulations of four different cosmological models, by identifying z ~ 3 
LBGs with the most massive halos in our simulations at that redshift. We are 
currently working on improving this analysis, both by comparing to the increasing 
amounts of new data, and by using semi-analytic models to explore different ways 
of identifying LBGs with halos in N-body simulations. 

HOW PROBABLE ARE THE SPIKES? 

Our methods are described fully in Wechsler et al. (W98) [1], and will be de- 
scribed only briefly here. The data is taken in 9'x9' fields; the redshift distribution 
is binned in bins of z — 0.4. In determining the probability of spikes, our statistics 
consider each pixel separately, where a pixel is one z=0.4 redshift bin by the size 
of an angular field: 9'xl8' (in the case of data from S98), or 9'x9', (in the case of 
the newer data from Adelberger et al. (A98) [3]). The number of galaxies in each 
pixel is divided by its selection function; we then consider the galaxy overdensity 
per pixel: 5 g = (N — N)/N. These data are then compared to N-body simulations 
by Gross et al. [8] of four cosmological models: SCDM (Q — 1, h — .5, a 8 = 0.67), 
CHDM (n = l,n cdm = 0.8, Q u = 0.2, h = .5,<r 8 = 0.72), OCDM (Q = 0.5, h = 
.6, a 8 = 0.66), ACDM (fi = 0.4, tt A = 0.6, h = .6, a 8 = 0.77), in a 75 h^Mpc box; 
57 million CDM particles (+ 113 million HDM particles in CHDM). The halos 
are identified as virialized regions at z ~ 3. To assign LBGs to halos we assume 
that one LBG resides in each massive halo, then choose a mass cutoff to match 
the observed number density of candidates and randomly select 40% of these (be- 
cause S98 find redshifts for ~ 40% of candidates). The mass cutoff ranges from 
3 x 10 11 for CHDM to 9 x 10 11 for SCDM. We then "observe" pixels the size of the 
observational determined pixels, and compared the statistic 5 g to that calculated 
from the data. The cumulative probability distribution of this galaxy overdensity 
statistic is plotted in Figure 1, both for the analysis done in W98, which compared 
only to the data from the one 9'xl8' field published by S98, and compared to the 
more complete data set from six 9'x9' fields published by A98. We find that the 
probability of observing a spike as large as the largest seen by S98 (corresponding 
to 5 g = 2.6) is about 37%, 31%, 27%, and 6% for CHDM, OCDM, ACDM, and 
SCDM, respectively. The reason for the small probability for SCDM seems to be 
due to the shape of the power spectrum [1], [3]. The distribution of the newer data 
set published in A98 looks relatively well fit by all of the models; a full analysis 
will be published elsewhere. In future work, we will also take into better account a 
more accurate selection function for the data. 

BIAS & CORRELATION FUNCTION 

The clustering of LBGs is found to be biased with respect to the underlying dark 
matter. The bias factor b, defined as b = 5 g /5 m , has an average value of about 2-5 
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FIGURE 1. (a) The cumulative probability distribution of 8 g for the S98 data as well as for the 
massive halo distribution in the four models considered here. The distribution reflects the relative 
excess of halos in each pixel. Spikes the size of the largest spike found by S98 are found with 
reasonable probability (~ 25 — 40%) in all models except SCDM, and with marginal probability 
in SCDM. (b) The cumulative probability distribution of 6 g for the newer data set of A98 as well 
as for the distribution of massive halos in the four cosmologies considered. The distribution seems 
reasonably well matched by all of the models. 



depending on the cosmological model. High-density regions are the most biased. 
We have measured the correlation function for the LBGs in our simulations. The 
best fit parameters are: r [h^ 1 Mpc] = 3.3, 5.1, 5.0, 7.3, and 7 = 1.7, 1.6, 1.6, 
1.5, for SCDM, CHDM, OCDM, and ACDM, respectively. A98 have calculated the 
correlation function using a counts-in-cells method, and find correlation lengths of 
~ 4 ± 1,5 ± 1, and 6 ± l(/i _1 Mpc) for Qm — 1 ; 0.2 open, and 0.3 flat, assuming 
that 7 = 1.8 (cf. [9]). 



ADDING SEMI-ANALYTIC MODELS 



In previous work, we identified LBGs by making the assumption that one object 
resides in each massive halo above some mass cutoff. Alternatively, we can use 
semi-analytic models [10], [11], [12] to predict the location of objects luminous 
enough to be observed as LBGs - they may reside in less massive halos, or there 
may be more than one per halo. We are currently exploring how different models of 
galaxy formation affect the clustering properties of the Lyman-break galaxies when 
identified in this way - which may help to distinguish between galaxy formation 
models and to determine the nature of the Lyman-break galaxies. 



CONCLUSIONS 



Large peaks in the observed redshift distribution of LBGs are common at z ~ 3 
in several cosmological models. Galaxy formation is biased at high-z in all cosmolo- 
gies we have considered, with a bias factor of ~ 2 — 5. The clustering properties 
will probably not distinguish between different cosmologies, without independent 
information about the bias. We find a similar correlation length but slightly shal- 
lower correlation function slope compared with that observed. We find that the 
bias and clustering properties are more strongly affected by the shape of the power 
spectrum than by the mass density; a model with a shallow power spectrum, like 
SCDM, seems to be somewhat less clustered than the data. 

A preliminary analysis of the data published by A98 shows that the original 
spike [2] is one of the highest in all surveyed fields, and the new data is in fairly 
good agreement with all four cosmological models. Major improvements to our 
prior analysis can be made both by comparing to the still-increasing amounts of 
new data, and by improving our method of assigning LBGs to dark-matter halos. 
By using semi-analytic models to test and improve upon our method of identifying 
LBGs within halos from N-body simulations, we may be better able to understand 
the clustering properties and nature of the Lyman-break galaxies, and may begin 
to distinguish between galaxy formation scenarios. 
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